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ABSTRACT 

Using the Cosmic Background Imager, a 13-element interferometer array operating in the 26-36 GHz 
frequency band, we have observed 40 deg 2 of sky in three pairs of fields, each ~ 145' x 165', using 
overlapping pointings (mosaicing). We present images and power spectra of the cosmic microwave 
background radiation in these mosaic fields. We remove ground radiation and other low- level contami- 
nating signals by differencing matched observations of the fields in each pair. The primary foreground 
contamination is due to point sources (radio galaxies and quasars). We have subtracted the strongest 
sources from the data using higher-resolution measurements, and we have projected out the response 
to other sources of known position in the power-spectrum analysis. The images show features on scales 
~ 6 / -15 / , corresponding to masses ~ 5-80 x 1O 14 M at the surface of last scattering, which are likely 
to be the seeds of clusters of galaxies. The power spectrum estimates have a resolution A£ « 200 and 
are consistent with earlier results in the multipole range i < 1000. The power spectrum is detected 
with high signal-to-noise ratio in the range 300 < £ < 1700. For 1700 < £ < 3000 the observations are 
consistent with the results from more sensitive CBI deep-field observations. The results agree with 
the extrapolation of cosmological models fitted to observations at lower i, and show the predicted 
drop at high £ (the "damping tail"). 

Subject headings: cosmic microwave background — cosmology: observations — techniques: inter fero- 
metric 



1. INTRODUCTION 

The Cosmic Background Imager (CBI) is a 13-element 
radio interferometer array designed to image the cosmic 
microwave background radiation (CMB) and measure its 
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angular power spectrum in the 26-36 GHz frequency 
band. The power of accurate measurements of the CMB 
power spectrum to constrain cosmological models and 
obtain precise estimates of critical cosmological parame- 
ters has been demonstrated by many theoretical and ob- 
servational studi es. The most recent e xperiments - the 
BOOMERANG (iNetterfield et, al.ll20?rl and MAXIMA 
(|Lee et al.ll2001j) balloon-borne bolomet ers and the DASI 
interferometer array at the South Pole (Halverson et al. 
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2002) - have measured the power spectrum at multipoles 
£ < 1000 (angular scales > 20'). The CBI has the poten- 
tial to extend these measurements to £ ~ 3500. 

This paper is the third in a series reporting results 
from the CBI. P reliminary results were presented by 
Pa din et, al.l (I200TI hereafte r Paper I). The accompanying 
paper ([Mason et all [2003. hereafter Paper II) presents 
our estimate of the power spectrum for £ < 3500 from ob- 
servations of three pairs of 45' (FWHM) deep fields made 
in our first observing season, 2000 January-December. 
The present paper (Paper III) presents complementary 
results from first-season observations of three pairs of mo- 
saic fields each about 145' x 165' (total « 40 deg 2 ), using 
the mosaicing method. These observations have higher 
resolution in £ than those of Paper II, and they have 
greater sensitivity at low £ owing to the reduced cosmic 
variance, but they are less sensitive at high £. Further 
observations made in 2001, which are currently being an- 
alyzed, will increase the sensitivity and improve the res- 
olution of our power spectrum estimate. The method we 
use for extracting power spectr um estimate s from inter- 
ferometry data is described bv lMvers et alJ (|2003l here- 
after Paper IV). 

This paper is organized as follows. In § we sum- 
marize the important properties of the CBI and intro- 
duce the mosaic technique. In § 03 we describe the ob- 
servations that are presented in this paper and present 
images of the three pairs of mosaic fields. In § ^ we 
describe the maximum-likelihood method for estimating 
the power spectrum from visibility measurements and 
present our power-spectrum estimates. We pay particu- 
lar attention to the contaminating effects of foreground 
point sources. Finally, in §[H|we discuss some of the im- 
plications of our results and summarize our conclusions. 
A full discussion of the implications for cosmology will 
be the subject of t wo further papers ([Sievers et al.l l2003. 
hereafter Paper V: lBond et al.ll2003l hereafter paper VI). 

2. THE COSMIC BACKGROUND IMAGER 

Theoretical models of the CMB predict its angular 
power spectrum, 

Ct = <k m | 2 ), (1) 

where a£ m are the coefficients in a spherical-harmonic ex- 
pansion of the CMB temperature distribution as a func- 
tion of direction a?, 



T(x) 



(2) 



£=0 1 



and 



2. 725 K is the mean CMB temperature 
1 Il999f) . The angle brackets indicate the 
expectation value (ensemble average). In this paper, as 
in most work, the quantity presented in the figures is the 
power per unit logarithmic interval in £, 

Ci = £(£+l)Ci/2ir, (3) 

scaled by Xq to put it in temperature units (/iK 2 ). 

Measurement of the CMB power spectrum with inter- 
ferometers has been discu s sed in several papers (e.g., 
Hobson. Lasenbv. fc Jonesl 11995 1: iHobson fc M aisinger 
20021 iMaisinger. Hobson. fc Lasenbvl 119971: iNeTEoOl: 
et al.lll999allbr ). and details of the method that we 
have used are presented in Paper IV. A single-baseline 
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Fig. 1. — Radial profile of the CBI primary beam in one of the ten 
frequency channels; data from all 78 baselines are superimposed. 
The observations (error bars) were made on 2000 Nov 12; Tau A 
was observed for about 20 s at each point on a 13 x 13 grid in 
azimuth and elevation. The grid points were separated by 7', so 
the grid extended to ±42'. Measurements within 45' of the central 
position were used to fit for the amplitude scale and pointing offset 
(3 parameters) for each of the 780 datasets (78 baselines x 10 
channels). The red curve shows the adopted profile, which was 
computed by taking the square of the Fourier transform of the 
aperture illumination pattern, assumed to be circularly symmetric. 
Taking the outer radius of the aperture (0.45 m) as r = 1, the inner 
part r < 0.172 is blocked by the secondary and is assumed to have 
zero illumination. The illumination is tapered from center to edge: 
we approximate this as a Gaussian exp[— (r/ro) 2 ]. The parameter 
7*0 = 0.683 was chosen to give the best fit to all ten channels. It 
corresponds to an edge taper of —18.6 dB and a beam FWHM of 
#fwhm — 45 .'2 x (31GHz/V). The 13 antennas have very similar 
beams, but they have relative pointing offsets of 1-2', up to 5 7 in 
the worst case. 



interferometer is sensitive to a range of multipoles £ ~ 
2ttu±A£/2 where u is the baseline length in wavelengths, 
and A£ is the FWHM of the visibility window function, 
which is proportional to the square of the Fourier trans- 
form of the primary beam (antenna power pattern). For 
a circular Gaussian primary beam of FWHM #fwhm rad, 

A^ = 4v / 21n2/^ FW HM. 

The CBI is a 13 element interferometer in which all 
78 antenna p airs are cross- corr elated (for a detailed de- 
scription, see llPadin et al. 2002). Its 26-36 GHz band is 
split into ten channels each 1 GHz wide, which are cor- 
related separately, giving a total of 780 complex visibil- 
ity measurements in each integration. The antennas are 
arranged with a common axis on a flat platform. The 
platform mount is steered in altitude and azimuth so 
that all the antennas track the same point on the celes- 
tial sphere; in addition, the platform is rotated about 
the axis to track parallactic angle, so that each baseline 
keeps a constant orientation relative to the field of view. 
The 78 baselines range in length from 1.0 m to about 
5.5 m, depending on the antenna configuration on the 
platform. During the observations reported here, we used 
several different configurations. The antennas respond to 
left circular polarization (LCP), although for part of the 
observations one antenna was configured for right circu- 
lar polarization (RCP). Data from the 12 cross-polarized 
baselines have not been used for this paper, but they 
will be used to place limits on CMB polarization (J. K. 
Cartwright et al., in preparation). 

The CBI Cassegrain antennas have a diameter of 
0.90 m and a measured primary beam width #fwhm = 
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45f2 x (31 GHz/z/) at frequency v, so that At « 300. The 
primary beam is quite close to a circular Gaussian, but 
we have adopted a more accurate model of the radial pro- 
file (see Fig.nj and used this model when making images 
and estimating the CMB power spectrum. 

In observations from a single pointing (as in Paper I 
and Paper II), the resolution in i of the power spec- 
trum is limited to ~ A£, which is insufficient to re- 
solve the expected structure in the spectrum. To im- 
prove the resolution in £, which is inversely proportional 
to the angular size of the imaged region, we make mo- 
saiced observations in which we map a larger area of sky 
using several closely-spaced pointings. This method is 
in widespread use for making images of extended re- 
gions with radio interferometers (e.g ., iCornwelJ 1988; 
Sault. St avelev-Smith. fc BrouwMl996T) . In the observa- 
tions reported here, we mapped three separate mosaic 
fields in this way, using 42 pointings for each in a rectan- 
gular grid of 7 rows separated by 20' in declination and 
6 columns separated by l m 20 s « 20' in right ascension. 
This allows us to improve the resolution in £ to A£ ^100. 

Although the CBI an tennas were design ed to have low 
sidelobes and crosstalk ()Padin et al.l200(J) . emission from 
the groun d contaminates th e data, especially on short 
baselines ()Padin et alJ 1200*3 Paper I; Paper II). The 
ground signal is stable on time-scales of many minutes, so 
if we observe two nearby fields under similar ground con- 
ditions, the difference of the visibilities of the two fields 
is unaffected by the ground. Differencing also eliminates 
any constant or slowly varying instrumental offsets. We 
observe a field (the lead field) for about 8 min and then 
switch to a reference field (trail field), at the same dec- 
lination but 8 min later in right ascension, for the next 
8 min, and form the difference of corresponding 8.4 s inte- 
grations. One such "scan" consists of up to 50 differenced 
integrations, the exact number depending on how much 
time is lost to slewing and calibration. The two fields are 
observed over the same range of azimuth and elevation, 
so they have nearly identical ground contributions. All 
the results presented in this paper are derived from the 
differenced visibilities. The images show the difference 
of intensity between a region of sky and one 8 min later 
in right ascension; and the differencing is included in the 
covariance matrices used for power-spectrum estimation. 

For each pointing incorporated in the mosaics we ob- 
tained approximately 16 such 8 min lead-trail scans, al- 
though the exact number varied from one pointing to an- 
other. Between scans, we rotated the platform to change 
its orientation relative to the hour circle, thus improving 
the sampling of the (u,v) plane 4 . This also reduces the 
effect of any residual ground contamination: the ground 
signal does not add coherently when the visibilities from 
different baselines measuring the same (u, v) point are 
combined, because the antennas have different far side- 
lobe responses. An example of the (u, v) plane sampling 
obtained for a single pointing is shown in Figure No 
attempt was made to obtain identical sampling for all 
the pointings. 



4 The vector u = (u, v) is the separation of a pair of antennas, 
measured in wavelengths (A), in a plane perpendicular to the di- 
rection of the center of the field of view, i.e., in the plane of the 
rotating antenna platform. 



3. OBSERVATIONS 

3.1. Summary of the dataset 

In this paper we present observations of three mosaic 
fields (identified as 02 h , 14 h , and 20 h ) separated by about 
6 h in right ascension at a declination of about —3? 5. The 
fields were chosen to have IRAS 100 /im emission less 
than 1 MJy sr _1 , low synchrotron emission, and no point 
sources brighte r than a few hundred mJy at 1.4 GHz 
((Condon et al.lll998l) . Details of the fields observed are 
given in Table 03 Two of the individual pointings com- 
prising our mosaics have been observed to much greater 
depth, and our data for these two pointings are a subset 
of the data analyzed in Paper II. 

Data acquisition, calibration, and editing were per- 
formed in the manner described in Paper II; we give 
only a summary here. Observations were made at eleva- 
tions > 42°, at night, and more than 60° from the moon. 
The amplitude scale was based on nightly observations 
of calibration sources (Jupiter, Saturn, Tau A [3C 144, 
the Crab Nebula], and Vir A [3C 274]); the primary 
calibrator was Jupiter, for which we assu med an effec- 
tive t emperature 5 of 152 ± 5 K at 32 GHz (M ason et alJ 
1999). We estimate that the overall calibration uncer- 
tainty is 5% rms, equivalent to 10% in CMB power. A 
small fraction of the data were discarded owing to instru- 
mental problems. Most of these problems were detected 
by real-time monitoring of the receivers; a few hardware 
problems in the correlators were indicated by unusually 
high correlation between the real and imaginary parts of 
the visibility. A few nights were affected by bad weather, 
and we deleted all data taken at times when atmospheric 
noise was visible on the short baselines. 

As mentioned above, all the data were taken in pairs of 
8 min scans on a lead and a trail field, separated by 8 min 
in RA. The field separation on the sky varied with decli- 
nation, but for our declinations, ~ — 3?5, it was very close 
to 2°. Individual 8.4 s integrations in the two scans were 
matched in hour angle and differenced, and unmatched 
integrations were discarded. We estimated the noise in 
each scan from the rms of the differenced integrations, 
and compared it with the expected rms noise level in 
either the real or imaginary visibility: 

a = ^L, (4) 
At]q v Avr 

where &b is Boltzmann's constant, T sys is the sys- 
tem temperature, A is the effective area of each 
antenna, tjq is the correlator efficiency, Av is the 
channel bandwidth, and r is the integration time 
((Thompson. Moran. fc Swensonll200lD . Under good con- 
ditions the measured noise, a « 4.7 Jy s 1 / 2 , is consistent 
with the estimated system temperature ~ 30 K. The rms 
noise in the differenced data should be \[2g. When the 
rms exceeded 2.6 times the expected rms, we discarded 
the entire scan-pair. This eliminated almost all of the 
data affected by the atmosphere. A final visibility esti- 
mate for each channel of each baseline in each orientation 
of the antenna platform was computed by a weighted 
average of the individual scan visibilities, and the un- 
certainty in this estimate was computed from the scan 

5 This is the excess brightness over the CMB, expressed as a 
temperature using the Ray leigh- Jeans approximation. 
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Fig. 2. — Typical (u,v) sampling obtained for a single pointing (C1444— 0230). For this observation only 12 antennas were used. Left: 
sampling for a single frequency channel. Right: sampling with 10 channels; a separate dot has been used for each channel. The different 
channels are sensitive to slightly different angular scales. 



rms's, taking into account a bias introduced by the fact 
that the scan rms's are themselves estimated from the 
data (see Paper II). 

The final edited and calibrated dataset is summarized 
in Table 03 which reports the total integration time (av- 
eraged over baselines) on each of the lead and trail fields 
comprising the mosaics. Owing to the vagaries of the 
weather and the telescope, some fields were observed to 
greater depths than others, and a few were missed al- 
together; this is reflected in the variation of sensitivity 
across each mosaic (see § E3> - 

3.2. Foreground Point Sources 

In the 26-36 GHz band, the dominant confusing fore- 
ground is the emission from discrete radio galaxies and 
quasars, which we refer to as "point sources" (they are 
virtually unresolved by the CBI). The contribution of 
point sources to the visibilities must be removed in order 
to obtain a reliable estimate of the CMB power spectrum. 
A random distribution of point sources has a power spec- 
trum Ci — constant, while for the CMB Ci decreases 
rapidly with increasing i, so discrete sources dominate 
at high I. 

To remove most of the point-source contamination 
in our data, we measured the flux densities of a 
large number of known point sources in the mosaic 
fields using a new dual-beam 31 GHz HEMT re- 
ceiver with a beamwidth of 81" on the 40-meter tele- 
scope at the Owens Valley Radio Observatory (OVRO). 
The observations followed standard methods (see, e.g., 
M vers. Readhead. fc Lawrencel |1993); details will be 
presented elsewhere (B. S. Mason et al., in prepara- 
tion). A total of 2225 sources brighter than 6 mJy at 
1.4 GHz in the NRAO VLA 1.4 GHz Sky Survey (NVSS; 
iCondon et aTIIT998h were observed, as described in Pa- 
per II. For each source detected above a 4a threshold, 
the expected visibility (assuming a point source at the 
NVSS position with the OVRO flux density attenuated 
by the CBI primary beam) was subtracted from the CBI 



visibility data. The survey is 90% complete for S31 gh z > 
16 mJy and 99% complete for 631GHz > 21 mJy. A to- 
tal of 70 sources were subtracted from the 02 h data, 63 
from the 14 h data, and 68 from the 20 h data. To es- 
timate the response in each CBI frequency channel, we 
used the two-point spectral index a (1.4-31 GHz) for 
steep-spectrum sources (a < —0.5 where S v ex v a ) and 
an average spectral index a = —0.23 for the remainder, 
which may be variable; but as we are only extrapolating 
±16% in frequency the results are not very sensitive to 
the choice of a. 

By making images before and after source subtraction 
(see § 13. 3|) we have verified that the CBI and OVRO 
measurements are consistent. The resulting images are 
dominated by the CMB. However, source subtraction is 
not sufficient for accurate power-spectrum estimation, 
because there may be small residuals and there will also 
be unmeasured sources that have not been accounted for. 
When estimating the CMB power spectrum, we have 
adopted a more powerful method than source subtrac- 
tion: we have projected out th e know n p oint sources 
f following lHalverson et"aIll2002D (see S Although 
we subtracted the sources measured at OVRO, we also 
projected them out, so errors in the OVRO flux-density 
measurements should not affect the power spectrum es- 
timates. 

3.3. Images 

The quantity of primary cosmological interest, the 
power spectrum, is best estimated directly from the vis- 
ibility data, as discussed in §0] below. However, we can 
also make images of the CMB (convolved with the instru- 
mental point-spread function) by Fourier-transforming 
the visibilities. The images provide a good check for 
the presence of non-gaussian features in the CMB or in- 
strumental errors in the data, such as calibration errors 
(which would show up as residuals after point-source sub- 
traction). 

We have used standard aperture-synthesis techniques 
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(Tavl or. Carilli. fc Perlevl 119991) to make images from 
each pointing. The images are formed from linear combi- 
nations of the measured visibilities and we have not done 
any deconvolution or "cleaning." The "dirty" image Id 
is the Fourier transform of the sampled visibilities, 



w k [V^ cos[27rn fc • x] - V k sm[2iru k • x}} 



Id(x) 

(5) 

and is the convolution of A(x)I(x) (the primary beam 
response times the sky brightness) with the point-spread 
function or dirty beam Bd(x): 



B D (x) 



^2 k W k COs[27Tlifc • X] 



(6) 



In practice, we resample the visibilities on a grid and use 
an FFT t o compute the im age and beam using standard 
software (Shepherd 1997). The sky coordinates x are 
the Fourier conjugates to the baseline components u and 
correspond to direction cosines relative to the pointing 
center. The weights w k are usually chosen to be the sta- 
tistical weights ("natural weighting"), w k = l/cr k where 
a k is the standard deviation of the real or imaginary part 
of the complex visibility, estimated as described above. 
With this weighting, the variance of the dirty image is 



s 2 = ([I D (x) - (Id(x))] 2 ) = 



2 ' 



(7) 



A "linear" mosaic image Im( x ) is formed from several 
pointings by shifting each dirty image to a common phase 
center, correcting it for its primary beam, and making a 
weighted sum at each pixel: 



Im(x) 



E p W p (x)I D Ax)/A p (x) 

EMx) 



(8) 



The weighted mean gives a maximum-likelihood estimate 
if the weights are equal to the inverse variances of the 
corrected images, i.e., W p (x) = A 2 (x)/s 2 , where s 2 is the 
variance of the dirty image made from the pth pointing. 
Thus 



Im(x) 



E p A 2 (x)/sl 

and the variance in the mosaic image is 

1 



E v A 2 (x)/sl 



(9) 



(10) 



(ICornwel]lH98a ISault. St avelev- Smith, fc Brouwl [l996). 
The variance varies across the image, and it is necessary 
to truncate the image where the variance becomes ex- 
cessive. 6 We have used the model for the primary beam 
described in the caption to Figure d set to zero for radii 
> 90 / . In making the mosaic images, and in the power 
spectrum analysis described below, we have made the 
approximation that the sky is flat over the ~ 3° image, 
which is equivalent to an error of < 5° of phase on the 
longest baselines. 

6 The full covariance matrix of the image pixels can be calculated 
in a similar way, but as we estimate the power spectrum from the 
visibility data directly we do not need this for our analysis. 



Figures EHHI show mosaic images made from the CBI 
data. We can make a variety of different images by se- 
lecting different subsets of the data and adjusting the 
weights w k . 

Figures EI an d [HI show mosaic images made from 
the entire dataset using natural weighting, and similar 
images made from the visibility data after subtraction of 
the point sources measured at OVRO. Figure shows the 
variation of noise level across these images, and indicates 
the pointing centers. The effectiveness of the OVRO 
source subtraction can be seen by comparison of the "be- 
fore" and "after" images. Most of the sources have been 
removed successfully, although there are of course resid- 
uals owing to measurement errors or source variability, 
and a few sources can be seen that were not detected at 
OVRO (these sources were projected out in the power 
spectrum analysis). 

To investigate the sources further, we have made 
higher-resolution images (not shown) using only base- 
lines longer than 250A, and searched for peaks exceed- 
ing 5<tm in images of signal-to-noise ratio. All these 
peaks are coincident within 2' with sources that have 
Si a GHz > 3.4 mJy in the NVSS catalog. Of the 42 
sources detected, 37 are identified with NVSS sources 
with 5i.4GHz > 6 mJy which were detected at OVRO. 
The remaining five objects are associated with NVSS 
sources that were not detected at OVRO; one of them, 
with Si. 4 GHz < 6 mJy, was not observed at OVRO. As 
we reported in Paper II, we find that the number N of 
sources greater than flux density S at 31 GHz is 

-1.0 



N(S) ^2.8±0.7deg" 



(11) 



10 mJy / 

for 5 < S < 50 mJy. We use this result in § 14.21 to 
estimate the contribution to our power spectra of sources 
below Si. 4 GHz > 3.4 mJy (stronger sources are treated 
individually) . 

Figures [7[ |H1 and [31 show lower-resolution images made 
by reducing the weight of long baselines with a Gaus- 
sian taper; this suppresses the high-frequency noise and 
emphasizes the CMB emission. These figures include im- 
ages made from the upper and lower halves of the CBI 
frequency band. There is good visual correspondence 
between the two halves of the band; but note that the 
spatial-frequency sampling is not the same in the two 
halves, so a quantitative comparison is difficult. The vi- 
sual agreement is echoed in the good agreement between 
the power spectr a obtained from the two halves of the 
band (see S 14.3.4)) . These low-resolution images are dom- 
inated by sky signal rather than noise. Note, however, 
that the images are differences of two sky patches, and 
they are also missing the lowest spatial frequencies, so 
a comparison with images obtained with another instru- 
ment may be difficult. By inspection of signal-to- noise 
ratio images (computed using Eas.l51and rrU|) made from 
subsets of the data, we have found that signals corre- 
sponding to t < 1000 are detected with high significance 
(> 3<j), and there are some significant detections of indi- 
vidual features at higher t. The detected features range 
in angular size from ~ 6' to ~ 15', corresponding to mass 
scales at the surface of last scattering of ~ 5 x 10 14 to 
8 x 10 15 Mq, so these features are likely to be the seeds 
that would evolve into clusters of galaxies by the present 
epoch. 
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Fig. 3. — Images of the 02 h mosaic. Left: raw data; right: after subtraction of sources measured at OVRO. The images show the 
difference of the emission in the lead and trail fields. The coordinates are J2000. The right ascension scale applies to the lead field; add 8 m 
to obtain the right ascension of objects in the trail field. The same gray-scale range has been used for both images, and it does not show 
the full brightness range of the discrete sources. Light (positive) spots are discrete sources in the lead field, while dark (negative) spots are 
discrete sources in the trail field. The brightest source in this image has a flux density of about 67 mjy at 31 GHz. These images were 
made from the entire dataset using natural weighting, and they have been corrected for the primary beam response as described in § 13.31 
They have a resolution (FWHM) of 5'.2-5'.5 (the resolution varies slightly across the image, depending on the (it, v) coverage obtained for 
each pointing). 
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Fig. 4. — Images of the 14 h mosaic. Left: raw data; right: after subtraction of sources measured at OVRO. For details, see caption to 
Fig. El 



4. POWER SPECTRUM 

4.1. Algorithm 

Our algorithm for the estimation of power spectra from 
mosaic visibility data is described in Paper IV. We model 
the power spectrum as flat in each of a set of Nb con- 
tiguous bands, i.e., 



(with £q = 0), and take the band-powers q B (B = 
i,...,7V#) as a set of unknown parameters to be de- 
termined by maximizing the likelihood (the probability 
of obtaining the measured visibilities, if the model were 
correct, for given values of the parameters). If the signal 
and noise obey Gaussian statistics, which we assume, the 
likelihood is given by the multivariate Gaussian distribu- 



te = l{l + l)C £ /2n = q B , t B -i <£<£b (12) 
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Fig. 5. — Image of the 20 h mosaic. Left: raw data; right: after subtraction of sources measured at OVRO. For details, see caption to 



Fig. El 






Fig. 6. — The variation of rms noise level om across each mosaic image (see equation 1101 . The pointing centers are marked with x. 
The outer parts of the images, where the noise level exceeds a specified threshold, have been blanked. In the 02 h and 20 h mosaics, the 
minimum noise level is 2.0 mjy/beam and the blanking threshold is 10 mJy/beam. The corresponding numbers are 1.7 and 8.3 mJy/beam 
for the 14 h mosaic. The sky areas covered within these thresholds are 7.45 (02 h ), 7.13 (14 h ), and 7.20 (20 h ) deg 2 in each of lead and trail. 



tion for complex variates, 
1 



exp 



(-V+CT 1 



V 



(13) 



7r n det C 

where V is a column vector containing the complex vis- 
ibility measurements, C = (VV^) (a function of the 
parameters qs) is the covariance matrix of the visibil- 
ities, and t denotes the Hermitian conjugate. (Although 
eq. ^Hl is expressed in terms of complex visibilities, it is 
easier in practice to treat the real and imaginary parts 
of the visibilities as a double- length real vector). 
The correlation matrix C can be written as 



C = C N + ^ te C| + g res C r 



+tovRoC OVRO + g NV ssC^-, (14) 

where C N is the noise correlation matrix, estimated from 
the data as described above, C% is the CMB signal cor- 
relation matrix for band B (independent of qs), C OVRO 



-iNVSS 



and C NVSS are constraint matrices (|Bond. Jaffe. fc Knoxl 
1998) representing the effects of foreground point sources 
of known position, and C res represents a residual con- 
tribution from faint sources of unknown position. We 
discuss the source terms further in § 14.21 The factors 
#res, #ovRO, an d #nvss could in principle be regarded 
as free parameters to be determined by maximum likeli- 
hood, but in practice they are not well determined by the 
data and we instead hold them fixed at a priori values. 

In a typical mosaic observation, the number of distinct 
visibility measurements is very large (~ 200 000 complex 
visibilities for a mosaic of 42 pointings with 10 frequency 
channels, 78 baselines, and several parallactic angles), 
which makes the covariance matrices (eq. |T3| ) imprac- 
tically large. However, neighboring points in the \u, v) 
plane are highly correlated and need not be treated com- 
pletely independently. To reduce the size of the matri- 
ces, we interpolate the measured visibility set V onto 
a smaller number of grid points in the (u,v) plane (so 
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that the quantities that enter the likelihood calculation 
are linear combinations of the measured visibilities), and 
make the corresponding transformation of the covariance 
matrix. By this means we reduce the dimension of the 
(real) matrices to < 5000. We ran tests with different 
grid spacings, using both real data and simulated data, 
to verify that using a finer grid would not significantly 
change the results (see Paper IV). We find the maximum 
li kelihood solution by the qua dratic relaxation technique 
of lBond. Jaffe. Knoxl JT998), which yields estimates of 
the band powers qs , with their covariances given by the 
inverse of the Fisher information matrix 

We also calculate the band-power window functions and 
the equivalent band-powers of the noise, known point 
sources, and residual point sources by the methods de- 
scribed in Paper IV . The band-power window function 
W B (£) (|Knoxlll999ft allows the expected band -power for 
a given model power spectrum to be estimated as a 
weighted mean: 

W=I >2^>. (16, 

4.2. Correction for Foreground Point Sources 

If we had accurate measurements of all the point 
sources in our frequency band, we could subtract their 
contributions directly from the visibilities, but as we do 
not we must adopt a statistical approach. A detailed de- 
scription of our method of dealing with point sources is 
presented in Paper II, and here we give only a summary. 

The strongest sources ("OVRO" sources) were mea- 
sured at OVRO and subtracted from the visibility data as 
described above. The subtraction was necessarily imper- 
fect, however, and the residuals of the subtracted sources 
make a contribution to the power spectrum that must 
be removed. We created a second list of "NVSS" sources 
within about 60' of any of the mosaic pointing centers 
that were not detected at OVRO but had 61.4 gh z > 
3.4 mJy, which corresponds to S31 GHz > 0.25 mJy for 
a typical spectral index of —0.84. This list contained 
960, 918, and 974 sources in the 02 h , 14 h , and 20 h mo- 
saic fields. Crude estimates of the flux densities for these 
sources were subtracted from the visibilities as part of the 
power-spectrum estimation. For each mosaic, we con- 
structed two constraint matrices, C OVRO and C NVSS , 
using the positions and estimated residual flux-density 
uncertainties of the sources in the two lists. In principle, 
if our error estimates were correct, the constraint ma- 
trices would fully account for the source contributions 
and should be included in the maximum likelihood anal- 
ysis with prefactors govRO = #nvss = 1- However, after 
some experimentation, we decided to err on the side of 
caution and use large prefactors, govRO = #nvss = 10 5 , 
to give very low weight to modes that are affected by 
these known sources (much larger factors cause the co- 
variance matrix to be ill-conditioned). This is equiva- 
lent to marginalizing over th e unknown flux densities, or 
"projecting out" the s ources (|Bond. Jaffe. fc Knoxl 1998: 
lHalverson et al.l2002 N >. Our analysis is thus insensitive to 



Table 1. Band-Powers and Uncertainties 
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errors in the assumed flux densities of the sources, but 
at the cost of some loss in sensitivity (see Paper IV). 

Sources that are not included in the known-source lists 
also contribute power. We estimate the visibility covari- 
ance arising from such sources and include it as the resid- 
ual source term C res . The radio source counts and spec- 
tra used to compute this term are described in Paper II. 
As in that paper, we estimate that, for an NVSS flux- 
density cutoff at 51.4 GHz > 3.4 mJy, the amplitude of 
the residual correction is C£ es = 0.08 ±0.04 Jy 2 sr _1 (see 
Paper IV). The dividing line between known sources and 
sources included in the residual term is somewhat arbi- 
trary, so long as the residual term is computed correctly 
for the chosen flux-density cutoff. We have conducted 
tests to verify that our results are insensitive to the pre- 
cise choice of cutoff. For a full discussion of the source 
projection and the residual correction, see Paper II. 

4.3. Results 

4.3.1. Joint Mosaic Power Spectrum 

The primary result of this paper is the power spectrum 
of the CMB in the three mosaics treated jointly. For this 
analysis we have estimated the power spectrum in bins of 
width A£ = 200, with two alternate locations of the bins. 
The "even" binning has £ B = 200 + 2005 (1 < B < 16), 
while the "odd" binning has £b — 100 + 2005 (1 < B < 
16); here £b is the upper limit of the bin, as in equa- 
tion (|T2j) . In both cases the first bin is wider and starts 
at £ = 0. While we included bins at higher i, we report 
the results only for £ < 3000: at higher £, the mosaic data 
have very little sensitivity. The two sets of bins are of 
course not independent. The results are given in Tabled 
which gives for each bin the band-power the rms un- 
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Fig. 7. — Central part of the 02 h mosaic, imaged with lower resolution, after subtraction of sources measured at OVRO. A Gaussian 
taper has been applied to the visibilities, falling to 0.1 at a baseline length of 500 wavelengths. Left: all frequency channels (26-36 GHz); 
beam FWHM w 6.9 / . Center: high-frequency channels (31-36 GHz); beam FWHM w 6.6'. Right: low-frequency channels (26-31 GHz); 
beam FWHM « 7.1'. 




Fig. 8. — Central part of the 14 h mosaic, imaged with lower resolution, after subtraction of sources measured at OVRO. A Gaussian 
taper has been applied to the visibilities, falling to 0.1 at a baseline length of 500 wavelengths. Left: all frequency channels (26-36 GHz); 
beam FWHM w 6.9 / . Center: high-frequency channels (31-36 GHz); beam FWHM w 6.6'. Right: low-frequency channels (26-31 GHz); 
beam FWHM « 7.1'. 




Fig. 9. — Central part of the 20 h mosaic, imaged with lower resolution, after subtraction of sources measured at OVRO. A Gaussian 
taper has been applied to the visibilities, falling to 0.1 at a baseline length of 500 wavelengths. Left: all frequency channels (26-36 GHz); 
beam FWHM w 6.9 / . Center: high-frequency channels (31-36 GHz); beam FWHM w 6.6'. Right: low-frequency channels (26-31 GHz); 
beam FWHM « 7.1'. 



certainty in qs from the Fisher matrix, and the centroid 
of the window function £ e ^. The results are also displayed 
in Figure El an d the window functions are shown in Fig- 
ure [O] 7 With A£ = 200, the adjacent bins are anticor- 
r elated by about 16%. We have also computed power- 

7 The window functions and inverse Fisher matrices are available 
on the CBI web page, |http : //w ww . astro . caltech . edu/~t j p/CBI/ 



spectrum estimates using narrower bins with A£ = 140, 
for which the anticorrelation of adjacent bins is about 
24%, again using overlapping "odd" and "even" bins. We 
have used all four binnings (At, = 140 odd and even, and 
A£ = 200 odd and even) for cosmological-parameter esti- 
mation (see Paper V), and all four give consistent results. 
The component band-powers (defined in Paper IV) for 
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Fig. 10. — Joint power spectrum estimates for the three CBI 
mosaics. Band-power estimates have been made for two alternate 
divisions of the t range into bins: "even" binning {green squares) 
and "odd" binning {blue circles). The error-bars show ±1<7 uncer- 
tainties from the inverse Fisher matrix. Two minimal inflation- 
based models are shown. Red: fit to CBI plus COBE DMR; 
Q tot = l.o, Q b h 2 = 0.0225, tt cdm h 2 = 0.12, tt A = 0.6, n s = 0.95, 
r c = 0.025, Cio = 786 fiK 2 . Black: joint fit to CBI, DMR, 
DASI, BOOMERANG-98, VSA and earlier data; Q tot = 1.0, 
Q b h 2 = 0.02, £l cdm h 2 = 0.14, Q A = 0.5, n s = 0.925, r c = 0, 
Cio = 887 /iK 2 . For details, see Paper V. 



1500 
I 



Fig. 12. — Joint po wer s pectrum estimates for the three CBI mo- 
saics (the same as Fig. 1101 . The stars represent the power spectrum 
of the instrumental noise correction, and the triangles represent the 
power spectrum of the residua l so urce correction. The black curve 
is the joint model from Figure Hoi 




Fig. 11. — Window functions for the joint power spectrum esti- 
mates for the three CBI mosaics. Top: "even" binning. Bottom: 
"odd" binning. The vertical grey lines show the bin boundaries. 
Each function is normalized: W B {l)/l = 1. 



instrumental thermal noise (C N ) and the residual source 
correction (C res ) are shown in Figure El This figure 
shows that the residual source correction is negligible at 
£ < 2000; the thermal noise, however, exceeds the signal 
for £ > 1300 and is the dominant effect at high £. Both 
these corrections increase approximately, but not exactly, 
as £ 2 . The thermal noise depends on the sampling in the 
(u, v) plane and better-sampled £ bins have lower noise; 
while the residual source correction has a non-thermal 
spectrum and the magnitude of its contribution in any 
bin depends on the spectral sensitivity in that bin — this 
varies from bin to bin as the CBI frequency channels do 
not have the same (u, v) sampling. 

Figure shows two theoretical spectra for minimal 
inflation-based models with different parameters. As we 
discuss in Paper V, we have evaluated the posterior prob- 
abilities of the CBI and other datasets over a grid of mod- 
els in a 7-dimensional parameter space, using a variety 
of prior probabilities based on Hubble constant, large- 



scale-structure, and supernova-la observations. The first 
model displayed in Figure is the model that maxi- 
mizes the posterior probability of the CBI mosaic (using 
the "odd" A£ = 140 binning) and COBE-DMR results, 
with the weak-h prior on the Hubble constant; within the 
grid of models, this is also the best fit if we restrict the 
search to fiat models. The second model is intended to 
represent a current "concordance" model: it is the best 
fit of the CBI, DMR, DASI, and BOOMERANG-98 data 
with fiat and weak-h priors. We find the same model if 
we also include LSS, SN, or HST-/i priors. The param- 
eters for the two models are given in the figure caption. 
It is remarkable that the two spectra are so similar. The 
CBI data together with DMR place strong constraints 
on the allowed region of parameter space. These con- 
straints are consistent with those obtained from earlier 
CMB observations, even though the CBI is sampling an 
£ range a factor of two larger than that spanned by the 
earlier experiments. For a full discussion of parameter 
estimation from the CBI results, see Paper V. 

We compare our results with the earlier results from 
the BOOMERANG, DASI, and MAXIMA experiments 
in Figure In the region of overlap (300 < £ < 1000) 
the agreement is very good. A detailed comparison will 
be made in Paper V. 

We have also computed the power spectrum of the 
three mosaics using the wider bins chosen for analyzing 
the CBI deep field data (Paper II). The deep and mo- 
saic results are compared in Figure EI We have used six 
bins: £ = 0-500, 500-800, 880-1445, 1445-2010, 2010- 
3500, and £ > 3500. The last bin contains virtually no 
data for the mosaics, and we have omitted it from the 
figure. The first bin (£ < 500)) is poorly constrained 
by the data and has a large uncertainty owing to sam- 
ple variance. Two of the three deep fields lie within the 
area covered by the mosaics, so the results are not en- 
tirely independent. But if we ignore this complication, 
we can compare the two sets of band-powers by a chi- 
squared test, using the band-power covariance matrices. 
For this test, we assume that the likelihood function is 
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Fig. 13. — Compa rison of the joint power spectrum esti mates from the three CBI mosaics (Fig. 1101 with the measurements from 
BOOMERANG INetterfield et al.ll200l) . DASI (Halverson et al. 2002), and MAXIMA (Lee et al. 2001). The rectangles indicate the 68% 
confidence intervals on band-power; for BOOMERANG, the solid rectangles indicate the 68% confidence interval for the statistical and 
sample variance errors, while the hatched rectangles shows the amount by which a zblcr e rror in the beamwidth (12. / 9 zb 1'4) would shift 
the estimates (all up or all down together). The black curve is the joint model from Figure ITuI 
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Fig. 14. — Comparison of the joint power spectrum estimates 
from the three CBI mosaics (blue circles) with those from the three 
deep fields of Paper II (red squares); both have been computed for 
the t b ins used in Paper II. The black curve is the joint model from 
Figure ED 



approximately Gaussian and compute 

x 2 = e e - ^) (*e + *e) - 1 - $) 

B B' 

(17) 

where q^) , q^) are the band-power estimates and 

Mg£,, M BB , are the inverse Fisher matrices for the two 
datasets. Omitting the first and last bins, we find 
X 2 = 5.57, with 5 degrees of freedom. If the two data 
sets were drawn from the same population, a larger value 
would be obtained in 35% of trials, so we conclude that 
the two data sets are consistent. At I > 2000, where 
the deep observations show a significant signal, the mo- 
saic observations are less sensitive than the deep and are 



consistent both with the deep result and with no signal. 
In the bin 2010 < t < 3500, we find q B = 148 ± 203 
/iK 2 in the mosaics and qs = 510 ± 171 /iK 2 in the deep 
fields. The thermal noise band-power in this bin is much 
larger for the mosaic data set than for the deep data set, 
so the deep results are less sensitive to systematic errors 
in the noise estimation and are thus more reliable than 
the mosaic results. 

4.3.2. Peaks and Dips in the Power Spectrum 

The power spectra with bin size A£ = 200 (Fig. ITU|) 
suggest the presence of peaks and dips, but, owing to 
the ant icorrelat ions between adjacent bins, their reality 
is difficult to assess in this presentation, or in similar 
plots of the A£ = 140 bins. To assess their significance, 
we have searched for extrema in the power spectrum fol- 
lo wing techniques app l ied to the BOOMERANG data 
bv Ide Bernardis et alJ (|2002fl . For each triplet of adja- 
cent bins [i — l,i,i+ 1] we model the local band-power 
profile as a three-parameter quadratic form 

q B = a£ B -2b£ B J rC, £ = i - 1, i, i + 1, (18) 

where l B , £ B are band-average values of £ 2 and £. In 
' terms of the fitted parameters a, 6, c, the peak location 
is £ p k = b/a, its amplitude is C p k = c — b 2 /a, and its 
curvature ft p k = a. We have assumed that the measured 
<2b's are Gaussian-distributed with a covariance F BB ,. 
In this Gaussian approximation for the likelihood C(q B ), 
the likelihood £(a, 6, c) of the quadratic parameters is 
also a Gaussian: the maximum values a m ,6 m ,c m are a 
direct transform of the data band-power averages in the 
three bins, and the curvature at the maximum, which de- 
scribes the uncertainty in these parameter estimates, is 
simply related to Fb b> • The maximum-likelihood values 
of ^ p k, Cpk and ft p k are determined by a m , 6 m , c m , but the 
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Fig. 15. — Results for the peak and dip detection algorithm ap- 
plied to the A£ = 140 data. Peaks (red) and dips (blue) found in 
the odd binning are indicated by diamonds with ellipses showing 
the la confidence contours; those found in the even binning are 
shown by triangles with dashed confidence contours. Only detec- 
tions with curvature above la are shown. For comparison, the filled 
and open green squares show the peak and dip positions and am- 
plitudes predicted from earlier CMB observations in the context of 
inflation- motivated models: they are from ensemble averages over 
the Q-database we use for parameter estimation in Paper V (for 
the weak-/i prior). In t his figure the black curve is the CBI+DMR 
model from Figure HH1 



errors in these transformed variables are non-Gaussian. 
We estimate the errors by computing the local curvature 
of the likelihood near maximum, by a Jacobian transfor- 
mation. We consider a peak or a dip in the spectrum 
to be detected if two conditions are met: (1) the posi- 
tion ^ p k is within the range of multipoles covered by the 
band triplet; (2) the best fit quadratic has a curvature 
a m which differs from zero by at least la. 

The results of applying this algorithm to the A£ = 140 
odd and even binned data for £ < 2000 are shown in 
Figure El There are more detections in the odd binning 
than in the even binning, which is an indication that 
even the A£ = 140 bin size is larger than we would like 
for effective peak/bin detection. In the odd binning, we 
detect four peaks, at £ ~ 550,800,1150,1500, and four 
dips, at £ - 400, 700, 1050, 1400. The significance levels 
of the peaks, expressed as |a m |/a, are 2.1, 2.6, 1.8, and 
1.7, and those of the dips are 2.3, 2.4, 1.5, and 2.3. The 
two peaks and two dips detected with the even binning 
show good agreement with those detected with the odd 
binning. Maximum-likelihood values for the subthresh- 
old even-binned detections (with curvature < la) are 
in good accord with the odd-binned detections. Thus we 
have tentative detections of the second, third, fourth, and 
fifth acoustic peaks. With this dataset, even A£ = 140 
is barely fine enough to resolve the peaks, but we should 
be able to do better when we include the second season 
of CBI data. 

A second approach t o peak / dip detection was also used 
bv Ide Bernardis et alJ (|2002l) : given a class of theoreti- 
cal models with a sequence of peaks and dips, the sta- 
tistical distribution of positions and amplitudes can be 
predicted by ensemble-averaging over the full probabil- 
ity, the multidimens ional likelihood. We have used the 
same (^-database as de Bernardis et alJ ((2002), which is 
also the one we have used for cosmological parameter 
estimation in Paper V. Figure IT51 shows the peaks and 
dips we "predict" from BOOMERANG, DASI, MAX- 
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Fig. 16. — Power spectrum estimates for the three CBI mosaics, 
treated independently. The upper panel shows the results for the 
odd binning, and the lower panel shows the results for the even 
binning. In both binnings, the lowest bin is poorly constrained by 
the data. In this bin, in particular, the sample variance is large. 
The sample variance contribution to the error estimate is a fraction 
of the fitted band-power, and is thus unrealistically low for the 02 h 
mosaic which has low fitted band-power. 



IMA, COBE DMR and 19 other experiments predating 
this CBI dataset. The errors on the positions and heights 
determined this way are relatively small, comparable to 
the size of the symbols plotted. Within this set of mini- 
mal inflation-based models, the positions and amplitudes 
of the higher--^ peaks are largely determined by the posi- 
tions and amplitudes of the first few. It can be seen that 
the values found using our model-independent quadratic 
peak/dip-finder are in excellent agreement with the pre- 
dictions for £ < 1000. At higher i, perhaps there is a 
shift of peak placement, but we caution that our peak- 
position error bars, being derived from a Fisher matrix 
determined at the maximum likelihood, are only approxi- 
mate. Adding CBI to the rest of the pre-CBI experiments 
gives peak positions and amplitudes in good accord with 
those shown here, and indeed so does using just DMR 
and the CBI data. 

4.3.3. Individual Mosaic Power Spectra 

To check the consistency of the three mosaic data sets, 
and to look for variations in the power spectrum with di- 
rection on the sky, we have computed the power spectra 
of the three mosaics separately. The results are shown for 
the two alternate binnings in Figure El the very noisy 
points for i > 2000 have been omitted. The chi-squared 
test (eq. [T7) ) shows that the 02 h and 14 h mosaics are 
consistent with each other, but the 20 h mosaic is dis- 
crepant at the 95% confidence level in the odd binning 
(see Table EJ}. Figure [HD shows that most of the discrep- 
ancy occurs in two adjacent bins at 700 < £ < 1100. It 
is unlikely that this reflects a real change in the CMB 
spectrum with direction, or that it could arise from, for 
example, foreground contamination in the 20 h mosaic, 
so we provisionally attribute the discrepancy to a chance 
statistical fluctuation. The three mosaics have galactic 
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Table 2. Chi-Squared Comparison of the Mosaic Power Spectra 



Mosaics 


Binning a 


X 2 (d.o.f.) b 


Significance Level 
(per cent) c 


02 h -14 h 


even 


7.8 (9) 


56 


02 h -14 h 


odd 


5.6 (9) 


78 


02 h -20 h 


even 


10.3 (9) 


33 


02 h -20 h 


odd 


17.4 (9) 


4 


14 h -20 h 


even 


7.6 (9) 


57 


14 h -20 h 


odd 


17.5 (9) 


4 



a The two alternate binnings are not independent. See § 14.3.11 
The first bin has been omitted. The £ range is 400-2200 in the 
even binning and 300-2100 in the odd binning. 

b X 2 computed using the inverse Fisher matrix, assuming Gaussian 
likelihood. 



c Probability of exceeding the observed % 2 if the two data sets are 
drawn from the same population. 
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Fig. 17. — Joint power spectrum estimates for the three CBI 
mosaics, split into the upper and lower halves of the CBI frequency 
band. The upper panel shows the results for the odd binning, and 
the lower panel shows the results for the even binning. 



latitudes -53° (02 h ), 48° (14 h ), and -27° (20 h ). The 
good agreement between the three mosaics suggests that 
the CMB power spectrum is not heavily contaminated 
by latitude-dependent diffuse emission from the Galaxy. 
We will be able to examine the question of field-to-field 
consistency more closely when we have analyzed data 
from the 2001 observing season, which extend the sky 
coverage by a factor of two. 

4.3.4. Subdivision by frequency 

A second consistency check is to compare the power 
spectra obtained from the different CBI frequency chan- 
nels. As in 8 13.31 we have divided the data into low- and 
high-frequency halves (26-31 and 31-36 GHz). The re- 
sults are shown in FigureEl If the signal is due primarily 
to the CMB, the two spectra should be consistent, but 
if there is a large contribution from a non-thermal fore- 
ground, such as synchrotron, free-free, or dust emission, 
they should be different. The two spectra are similar, 
but it is difficult to make a quantitative comparison ow- 



ing to the strong correlation between the two frequency 
bands. The error estimates obtained from the Fisher ma- 
trix include the contribution of cosmic variance as well as 
the measurement noise, so they overestimate the band- 
to-band variance. 

5. CONCLUSIONS 

The results presented in this paper demonstrate the 
effectiveness of interferometric mosaicing both for imag- 
ing the CMB and for measuring its power spectrum. The 
CBI images show for the first time structures in the CMB 
on mass-scales corresponding to clusters of galaxies, and 
the CMB power spectrum has been extended by more 
than a factor of two in multipole number £. Although 
the resolution in £ is limited (we will obtain better reso- 
lution when we analyze the wider-field observations made 
in 2001), we have been able to detect the second and third 
acoustic peaks in the spectrum, and, for the first time, 
the fourth and possibly the fifth. The deeper observa- 
tions reported in Paper II extend the spectrum even fur- 
ther, well into the damping tail region where secondary 
anisotropics become important (see Paper VI). 

Ground-based observations of the CMB in the 1-cm 
wave-band, where long integrations can be obtained, are 
competitive with balloon-based observations at shorter 
wavelengths. The principal obstacle to observations in 
this band, particularly at high i, is the emission from 
foreground point sources. We have shown that it is pos- 
sible to correct the observations for this contamination 
with high accuracy, at a cost in sensitivity, but fore- 
ground sources remain the largest source of uncertainty 
in the power spectrum. To improve on our result, we will 
need sensitive high-resolution surveys of the foreground 
sources at 31 GHz, such as will soon be possible with the 
NRAO Green Bank Telescope. 

It should be clear from Figures ITOl and ESI that the 
CBI results are consistent with earlier observations in 
the region £ < 1000. What is more remarkable is that at 
higher i, a region that has not been probed before, the 
results are consistent with extrapolations of the power 
spectrum based on simple inflation- motivated models. 
We show in Paper V that the major cosmological pa- 
rameters (the total density parameter, the density pa- 
rameters for baryonic and non-baryonic matter, the pri- 
mordial density perturbation spectral index, the Hub- 
ble constant, the cosmological constant, and the optical 
depth to last scattering) are well constrained by the CBI 
observations, even when only the region of the spectrum 
beyond the first two peaks is considered. This provides 
further strong support for cosmological models domi- 
nated by cold dark matter and dark energy, and with 
a scale-invariant spectrum of primordial density fluctua- 
tions up to £ > 2000. The corresponding angular scales 
and masses are ~ 6' and ~ 5 x 10 14 M Q , the scale of 
galaxy clusters. This provides a firm foundation for the- 
ories of galaxy formation. 

A second season of CBI mosaic observations were ob- 
tained in 2001, and are currently being analyzed. These 
observations double the size of each of the three mosaic 
fields, and will enable the power spectrum to be deter- 
mined with improved sensitivity and ^-resolution. We 
are currently reconfiguring the instrument to maximize 
its sensitivity to polarization, with the goal of detecting 
and measuring the power spectrum of the polarized com- 
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ponent of the CMB, which is another powerful test of the 
cosmological models. 
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MOSAIC OBSERVATIONS WITH THE COSMIC BACKGROUND IMAGER 



Table 3. Summary of Observations 



Field name a 




Date 


it. A. 

(h, m, s) 


ueci. 


Integration 
time (s) 


C0242-0230 


2000 


Jul 16, Oct 20 


02 42 00 


-02 30 


4560 


C0242-0250 


2000 


Oct 01 


02 42 00 


-02 50 


7454 


C0242-0310 


2000 


Jul 31, Oct 21 


02 42 00 


-03 10 


4394 


C0242-0330 


2000 


Oct 04, Oct 06 


02 42 00 


-03 30 


14388 


C0242-0350 


2000 


Aug 03, Oct 21 


02 42 00 


-03 50 


5090 


C0242-0410 


2000 


Oct 18 


02 42 00 


-0410 


6170 


C0242-0430 


2000 


Aug 06, Oct 22 


02 42 00 


-04 30 


4954 


C0243-0230 


2000 


Sep 09, Oct 22, Oct 25, Oct 26 


02 43 20 


-02 30 


10220 


C0243-0250 


2000 


Aug 09, Oct 22 


02 43 20 


-02 50 


5048 


C0243-0310 


2000 


Sep 22 


02 43 20 


-0310 


6314 


C0243-0330 


2000 


Aug 29 


02 43 20 


-03 30 


5562 


C0243-0350 


2000 


Sep 25 


02 43 20 


-03 50 


5562 


C0243-0410 


2000 


Sep 02 


02 43 20 


-0410 


3206 


C0243-0430 


2000 


Sep 28 


02 43 20 


-04 30 


7462 


C0244-0230 


2000 


Oct 20 


02 44 40 


-02 30 


2052 


C0244-0250 


2000 


Oct 02 


02 44 40 


-02 50 


7576 


C0244-0310 


2000 


Aug 01, Oct 21 


02 44 40 


-0310 


4980 


C0244-0330 


2000 


Oct 07 


02 44 40 


-03 30 


6990 


C0244-0350 


2000 


Aug 04, Oct 21 


02 44 40 


-03 50 


5100 


C0244-0410 


2000 


Oct 19 


02 44 40 


-0410 


3706 


C0244-0430 


2000 


Aug 07, Oct 22 


02 44 40 


-04 30 


5052 


C0246-0230 


2000 


Sep 10 


02 46 00 


-02 30 


6812 


C0246-0250 


2000 


Aug 11, Aug 12 


02 46 00 


-02 50 


3480 


C0246-0310 


2000 


Sep 23 


02 46 00 


-0310 


6314 


C0246-0330 


2000 


Aug 30, Oct 24 


02 46 00 


-03 30 


9384 


C0246-0350 


2000 


Sep 26 


02 46 00 


-03 50 


6908 


C0246-0410 


2000 


Sep 07, Oct 25, Oct 26 


02 46 00 


-0410 


13380 


C0246-0430 


2000 


Sep 29 


02 46 00 


-04 30 


6766 


C0247-0230 


2000 


Jul 29, Oct 20 


02 47 20 


-02 30 


4854 


C0247-0250 


2000 


Oct 03 


02 47 20 


-02 50 


7574 


C0247-0310 


2000 


Aug 02, Oct 21 


02 47 20 


-03 10 


5060 


C0247-0330 


2000 


Oct 08 


02 47 20 


-03 30 


7380 


C0247-0350 


2000 


Aug 05, Oct 22 


02 47 20 


-03 50 


4910 


C0247-0410 


2000 


Oct 09 


02 47 20 


-0410 


2066 


C0247-0430 


2000 


Aug 08, Oct 22 


02 47 20 


-04 30 


5048 


C0248-0230 


2000 


Sep 11, Oct 26 


02 48 40 


-02 30 


8356 


C0248-0250 


2000 


Aug 15, Oct 23 


02 48 40 


-02 50 


9296 


C0248-0310 


2000 


Sep 24 


02 48 40 


-0310 


4404 


C0248-0330 


2000 


Sep 01 


02 48 40 


-03 30 


5134 


C0248-0350 


2000 


Sep 27, Oct 25 


02 48 40 


-03 50 


5826 


C0248-0410 






02 48 40 


-0410 




C0248-0430 


2000 


Sep 30 


02 48 40 


-04 30 


7170 


C1442-0230 


2000 


May 05 


14 42 00 


-02 30 


7288 


C1442-0250 


2000 


Jul 19, Aug 23 


14 42 00 


-02 50 


4054 


C1442-0310 


2000 


Apr 04, Apr 05, Apr 27 


14 42 00 


-03 10 


26048 


C1442-0330 


2000 


Jul 26, Aug 17 


14 42 00 


-03 30 


4560 


C1442-0350 b 


2000 


Mar 17, Apr 28 


14 42 00 


-03 50 


14610 


C1442-0410 


2000 


Aug 01, Aug 12, Aug 16 


14 42 00 


-0410 


6766 


C1442-0430 


2000 


May 01 


14 42 00 


-04 30 


7492 


C1443-0230 


2000 


Jun 24, Aug 19 


14 43 20 


-02 30 


6152 


C1443-0250 


2000 


May 23 


14 43 20 


-02 50 


5180 


C1443-0310 


2000 


Jun 27 


14 43 20 


-0310 


6814 


C1443-0330 


2000 


May 26 


14 43 20 


-03 30 


5210 


C1443-0350 


2000 


Jul 02 


14 43 20 


-03 50 


6402 


C1443-0410 


2000 


May 29 


14 43 20 


-0410 


7352 


C1443-0430 


2000 


Jul 15 


14 43 20 


-04 30 


5582 


C1444-0230 


2000 


May 11 


14 44 40 


-02 30 


8396 


C1444-0250 


2000 


Jul 20, Aug 18 


14 44 40 


-02 50 


5878 


C1444-0310 


2000 


Apr 9, Apr 10, Apr 26 


14 44 40 


-0310 


13920 


C1444-0330 


2000 


Jul 27, Aug 22 


14 44 40 


-03 30 


5134 


C1444-0350 


2000 


Apr 07, Apr 08, Apr 25 


14 44 40 


-03 50 


13704 


C1444-0410 


2000 


Jun 20, Jun 21, Jun 22 


14 44 40 


-0410 


20112 


C1444-0430 


2000 


May 02 


14 44 40 


-04 30 


6326 


C1446-0230 


2000 


Jun 25 


14 46 00 


-02 30 


5066 


C1446-0250 


2000 


May 25 


14 46 00 


-02 50 


5266 


C1446-0310 






14 46 00 


-0310 




C1446-0330 


2000 


May 28 


14 46 00 


-03 30 


4904 


C1446-0350 


2000 


Jul 03 


14 46 00 


-03 50 


6808 


C1446-0410 


2000 


May 30 


14 46 00 


-0410 


7366 


C1446-0430 


2000 


Jul 16 


14 46 00 


-04 30 


4826 


C1447-0230 


2000 


May 22 


14 47 20 


-02 30 


7254 


C1447-0250 


2000 


Jul 24, Jul 25 


14 47 20 


-02 50 


4002 


C1447-0310 


2000 


Apr 11, Apr 12, Apr 29 


14 47 20 


-0310 


22004 


C1447-0330 


2000 


Jul 29 


14 47 20 


-03 30 


4754 


C1447-0350 


2000 


Apr 30, Jul 23 


14 47 20 


-03 50 


8212 


C1447-0410 


2000 


Aug 15, Aug 21 


14 47 20 


-0410 


5166 
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Table 3. Summary of Observations — Continued 



Field name a 




Date 


IX. A. 

(h, m, s) 


ueci. 

/o /\ 

{ > ) 


Integration 
time (s) 


C1447-0430 


2000 


May 03 


14 47 20 


-04 30 


4606 


C1448-0230 


2000 


Jun 26 


14 48 40 


-02 30 


7088 


C1448-0250 


2000 


May 31 


14 48 40 


-02 50 


8490 


C1448-0310 


2000 


Jul 01, Aug 20 


14 48 40 


-0310 


5048 


C1448-0330 


2000 


Jun 04 


14 48 40 


-03 30 


6308 


C1448-0350 


2000 


Jul 04 


14 48 40 


-03 50 


818 


C1448-0410 


2000 


Jun 07 


14 48 40 


-0410 


8058 


C1448-0430 


2000 


Jul 17 


14 48 40 


-04 30 


5222 


C2042-0230 


2000 


Jun 04 


20 42 00 


-02 30 


2844 


C2042-0250 


2000 


Jul 28 


20 42 00 


-02 50 


5556 


C2042-0310 






20 42 00 


-0310 




C2042-0330 






20 42 00 


-03 30 




C2042-0350 






20 42 00 


-03 50 




C2042-0410 






20 42 00 


-0410 




C2042-0430 


2000 


May 05, Jul 02 


20 42 00 


-04 30 


9392 


C2043-0230 


2000 


May 31 


20 43 20 


-02 30 


5930 


C2043-0250 


2000 


Jun 08 


20 43 20 


-02 50 


7520 


C2043-0310 


2000 


Jul 05 


20 43 20 


-0310 


8024 


C2043-0330 


2000 


Jun 12 


20 43 20 


-03 30 


6992 


C2043-0350 






20 43 20 


-03 50 




C2043-0410 


2000 


Jun 26 


20 43 20 


-0410 


7128 


C2043-0430 


2000 


Jul 25 


20 43 20 


-04 30 


4356 


C2044-0230 


2000 


May 30 


20 44 40 


-02 30 


7536 


C2044-0250 






20 44 40 


-02 50 




C2044-0310 


2000 


May 01 


20 44 40 


-0310 


4536 


C2044-0330 






20 44 40 


-03 30 




C2044-0350 


2000 


May 02 


20 44 40 


-03 50 


2150 


C2044-0410 






20 44 40 


-0410 




C2044-0430 


2000 


May 11 


20 44 40 


-04 30 


4538 


C2046-0230 


2000 


Jul 03 


20 46 00 


-02 30 


8724 


C2046-0250 


2000 


Jun 10 


20 46 00 


-02 50 


6298 


C2046-0310 


2000 


Jul 06 


20 46 00 


-0310 


8722 


C2046-0330 


2000 


Jun 13 


20 46 00 


-03 30 


7684 


C2046-0350 


2000 


Jul 23 


20 46 00 


-03 50 


8302 


C2046-0410 


2000 


Jun 27 


20 46 00 


-0410 


8564 


C2046-0430 


2000 


Jul 26 


20 46 00 


-04 30 


8260 


C2047-0230 


2000 


Jun 07 


20 47 20 


-02 30 


7596 


C2047-0250 






20 47 20 


-02 50 




C2047-0310 


2000 


May 04 


20 47 20 


-0310 


3258 


C2047-0330 






20 47 20 


-03 30 




C2047-0350 


2000 


May 03 


20 47 20 


-03 50 


2764 


C2047-0410 






20 47 20 


-0410 




C2047-0430 


2000 


May 29 


20 47 20 


-04 30 


7132 


C2048-0230 


2000 


Jul 04 


20 48 40 


-02 30 


8312 


C2048-0250 


2000 


Jun 11 


20 48 40 


-02 50 


7326 


C2048-0310 


2000 


Jul 07 


20 48 40 


-0310 


8652 


C2048-0330 b 


2000 


Aug 02, Aug 23 


20 48 40 


-03 30 


15316 


C2048-0350 


2000 


Jul 24 


20 48 40 


-03 50 


8318 


C2048-0410 


2000 


Jul 01 


20 48 40 


-0410 


6660 


C2048-0430 


2000 


Jul 27 


20 48 40 


-04 30 


7630 



a The pointing center of the lead field is given; each is accompanied by a trail field 8 min later in R.A. Coordinates are J2000. 
b Data from this pointing were included in the deep dataset (Paper II). 



